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STRESS AND SHAPE ANALYSIS OF A PARAGLIDER WING 1 
By Robert W. Fralich^ 


SUMMARY 


^2 (Z $ Z-- 


A combined aerodynamic -structural analysis is made which is based on the 
assumption that the sail is flexible and has freedom to take the shape which the 
aerodynamic pressure and the internal stresses dictate. Numerical results were 
obtained for Newtonian impact aerodynamic theory and were compared with published 

results obtained for a rigid idealization of the paraglider wing. It was found 

•< 

that the assumed rigid idealization did not approximate the shape of a flexible 
wing well and led to significant errors in the lift and drag forces and the 
lift-to-drag ratio. The new calculations provide a basis for design of para- 
gliders for hypersonic flight. 

INTRODUCTION 

The paraglider wing is under active development as a controllable lifting 
device for landing space vehicles. It has also been suggested for other appli- 
cations such as recovering rocket boosters and effecting reentry into the earth’s 

atmosphere. A typical configuration is shown in figure 1 . ^See also [l]^ 

It consists of leading-edge booms and a keel boom joined together at the nose 

-klhis paper is based on a dissertation submitted in partial fulfillment of 
the requirements of the degree of Doctor of Philosophy at Virginia Polytechnic 
Institute, Blacksburg, Virginia, June 1963 . 

^Aero space Engineer, NASA Langley Research Center, Langley Station, 

Hampton, Virginia. 

^Numbers in brackets designate references at end of paper. 
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and a flexible sail whose surface carries the aerodynamic pressure loading. The 
payload is usually suspended beneath the wing by shroud lines. The distance 
between booms is often constrained by use of spreader bars or by means of a 
spring mechanism at the nose. Control is afforded by adjusting the length of 
the shroud lines and thus positioning the payload with respect to the wing. 

Various wind-tunnel and free-flight investigations have been made using 
either flexible or rigid idealizations of paraglider wings in order to find 
lift and drag characteristics or stability and control characteristics, (see, 
for example, [l-8j . ) None of the investigations on flexible sails yielded the 
deflected shape, the pressure distribution on the sail, or the stresses in the 
sail. Two investigations on rigid idealizations did yield pressure dis- 

tributions. However, these distributions apply to the fixed shape of the 
idealized wing and only represent the distribution on an actual flexible sail 
if the rigid idealization corresponds to the actual deflected shape of the 
flexible wing. 

In this paper a theoretical investigation is made which consists of a com- 
bined structural and aerodynamic analysis. In this analysis, the sail is 
assumed to be an inextensional, flexible membrane which has complete freedom to 
take the shape which the aerodynamic pressures and the internal stresses dic- 
tate. The assumptions are also made that the booms are rigid and straight and 
that they have small enough cross sections so as not to affect the aerodynamics. 
The booms are maintained a fixed distance apart by spreader bars, and the 
dihedral of the leading-edge booms is fixed with respect to the keel boom. 

The aerodynamic theory used is Newtonian impact theory, which has some- 
times been used to express the aerodynamic pres sure -shape relationship for the 
hypersonic speed range. (see, for example, jj9,lo].) This choice of aerodynamic 
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theory is governed by two factors. First, Newtonian theory leads to a simpli- 
fied analysis which shows the applicability of the present approach in a simple 
manner. Second, this aerodynamic theory has been applied in [ 7 ] to a rigid 
idealization of a paraglider wing, and thus serves as a basis of comparison for 
the numerical results. Numerical results are obtained for a flexible configu- 
ration that corresponds to this rigid idealization in the sense that it has the 
same surface planform and the same keel and leading-edge locations. The vari- 
ations of lift and drag coefficients with angle of attack are compared with 
those for the rigid idealization and the deflected shape, pressure distribution, 
and stress resultants in the sail are determined. Also, numerical results are 
presented to show effects of variation in dihedral angle (raising or lowering 
of the leading-edge booms) , 


GOVERNING EQUATIONS 
Geometry of the Sail 

In this analysis the sail is assumed to be constructed from a flat, inex- 
tensional membrane which may take some rather general shape in its loaded 
equilibrium condition. By virtue of its inextensibility, it is possible to 
establish appropriate coordinates in the surface by considering the sail in 
its flat (unloaded) condition. (See fig. 2(a).) Points on the sail are thus 
designated by the coordinates x and 0. Since the paraglider wing is sym- 
metric about the x-^ axis (fig. l) , and since only symmetric deformations and 
loadings will be considered herein, only that portion of the wing in the first 
quadrant need be considered. The keel boom, of length Z^-, is located at 
0=0 and the leading-edge boom, of length Z^, at 0 = 0^. It is also assumed 
that the trailing edge of the sail is straight; points on the trailing edge are 
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expressed by the equation 


= 




sin 0 + A cos 0 


(i) 


where 


sin 0 t 

A = 

l K 

— - cos 0^ 


( 2 ) 


and a subscript T has been added to indicate that the values of x from this 
equation are values at the trailing edge. 

For the x,0 polar coordinate system, the first fundamental form is given 

by 


ds 2 = dx 2 + x 2 d0 2 


( 3 ) 


Since the booms are straight and rigid and lie on x-coordinate lines and since 
deformations of the sail are assumed inextensional, the x,0 coordinate curves 
are lines of principal curvature of the loaded surface for which 


R i = “ 


Ro = — 








are the radii of principal curvature. This is a surface of zero Gaussian 
curvature which has the second fundamental form gd9^. Note that the sign con- 
vention on g differs from that usually employed in references in differential 
geometry of surfaces in that, for positive curvature of the surface, the posi- 
tive direction of the normal is outward. The only one of the Codazzi equations 
of compatibility (see [ll] ) not identically satisfied becomes 
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_ g 

dx x 


Integration of equation (5) yields 


( 5 ) 


g = 

R(e) 

where R(0) is an arbitrary function of 9 alone, 
ture R 2 (eq. (4)) is 


(6) 


Then the radius of curva- 


R 2 = xR( 9) 


(7) 


Equilibrium Equations 

Equilibrium of the loaded sail is governed by the equilibrium equations 
of membrane shells of zero Gaussian curvature ^see, for example, [l2 J)- In 
terms of the x,0 coordinates, these are: 


*“x,x + “xe,e + "x - N e - 0 ( 8 > 

^xe.x + N 0,e + ^xe = 0 <9) 


N 0 = R2X = x(RX) (10) 

for the case where only normal pressure forces X act on the surface of the 
sail. Here N x and N 0 are the normal stress resultants in the x and 9 
directions, respectively, and N x q is the shearing stress resultant. The 
comma followed by a subscript denotes partial differentiation with respect to 
this subscript. 
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Boundary Conditions 

It is appropriate to consider force boundary conditions along a general 
boundary contour C of the loaded surface. Thus, one may prescribe 


P L " h\ + 2l l Z 2 N x0 + hX 

P S = 2 l s l N x + ( 2 1 S 2 + 2 2 s l) N x9 + 2 2 s 2 N 9 


(ID 


where Pg and Pg are applied boundary forces per unit length in the surface 
normal and tangent to C, respectively. Here 7 ^ and 7g are l* 16 components 
of the unit outward surface normal L to C and s^ and S2 are the compo- 
nents of the unit tangent S (in the direction of increasing S) . Since L 
and S are orthogonal vectors 


s-^ 7 -| t S 2 ^2 — 0 

Sj_ 2 + 82^ = 1 
l ± 2 + * 2 2 = 1 


( 12 ) 

(13) 

(14) 


From figure 3 it can be seen that 



s 2 





(15) 


For a given boundary contour C, and s 2 can be found from equa- 

tions (13) and (15) and 7-^ and 72 from equations (12) and (l4). Then the 
force boundary conditions for the given contour are determined from 
equations (ll) • 
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Trailing-edge boundary conditions .- At the trailing edge the boundary con- 
tour is expressed by equation (l). From equations (l) and (15) the following is 
obtained : 

(sin 0 + A cos + (cos 9 - A sin 0)s 2 = 0 ( l6) 

Simultaneous solution of equations (13) and (l6) together with the solution of 
equations (12) and (l4) yields 

s n = -Ip = - = (cos 6 - A sin 9) (17) 

1 ATia 


** 


1 ~ /T7 


=(sin 0 + A cos 0) 


(18) 


Now application of equation ( 11) along a stress-free trailing edge yields 


and 


( p l) t = — ^ (sin 0 + A cos 6) 2 (n x ) t 
1 "J" A 


L 


+ 2(cos 0 - A sin 0)(sin 9 + A cos 9) ( N *e) T 
+ (cos 0 - A sin = ® 


(19) 


( P s) T = 1 — j( sin 9 + A cos 0)( cos e - A sin e)|(N 0 ) T - (Nx)t ] 

1 "J* A 

+ [(sin 0 + A cos 0) 2 - (cos 0 - A sin 0 )^J(n xQ ) t J» = 0 (20) 

Finally, equations (19) and (20) may be recast in convenient form by solving for 
(V) T and and using equations (l) and (10) 



( 21 ) 


and 



(RX) t 


2 

(RX) t 


( 22 ) 


Nose boundary conditions .- The boundary conditions at the nose of the sail 
are obtained by considering the limiting case of a stress-free boundary that 
shrinks to a point. Hence the stress resultants N x , Ng, and N x g must 
remain finite at this point . 


Keel and leading-edge boundary conditions .- The boundary conditions at the 
keel and leading edge are given by specifying the relative positions of the keel 
and leading edge. Details of the specification of these boundary conditions 
are reserved to a subsequent section. 


ANALYSIS 


Solution of Equilibrium Equations 

Substitution of the expression for Ng from equation (10) into equa- 
tions (8) and (9) yields 


(xN x ) + N x6;6 = x(RX) 


(23) 


(x 2 Nxe ) = -x 2 (RX) 


,9 


(24) 


Integration of equation (24) gives 


N. 


X0 


£ 1 


0 


(25) 
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where the boundary condition for finiteness of the stress resultant at the nose 
of the sail has been satisfied. Substitution of N x g from equation (25) into 
equation (25) and integration then yields 



, rt 

N x = - / 
x J 0 

§(RX) + -i- / A 2 (RX) dA 

i 2 J o ’ QQ 


(26) 


where again, the boundary condition at x = 0 has been imposed. Thus, the 
stress resultants Nq, N x q, and N x are readily expressed in terms of the 
aerodynamic pressure X(x,0) and the single parameter R(9). 

Now the boundary conditions (21) and (22) for stress resultants at the 
trailing edge of the sail are seen to require: 




e d ^ 



(27) 



*T J 0 



09 


dA 


as 


xpydo / 


(KX) r 


(28) 


t should be noted that both equations (27) and (28) must be satisfied; how- 
ver, only the one aribtrary function R(0) is available for satisfying these 
wo conditions. Difficulties of this nature are commonly found when dealing with 
lembrane shells of zero Gaussian curvature (see, for example, [12] ) and as a 
•esult membrane stress states can, in general, be obtained only in special 
cases. For the straight trailing-edge boundary conditions considered in this 
problem the fortunate situation arises that, for any aerodynamic pressures of 
the type 

00 

Y, a n xnz n( 0 ) 

n— 1 


X(x,0) = 


(29) 
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a single function R(e) can be found which satisfies both equations (27) 
and (28) and thus leads to true membrane stress states. Moreover, for aero- 
dynamic pressures of this type, a straight trailing edge is necessary for 
obtaining a membrane stress state that satisfies both trailing-edge boundary 
conditions. 

The parameter R( 8) expresses the shape of the loaded surface of the sail 
and the aerodynamic pressure X(x,0) is dependent on this shape. It is con- 
venient for deriving a relationship between aerodynamic loading and the shape 
to express the shape in terms of angles p(0) and 5(0) measured from the keel 
location rather than by R(0). (See fig. 2(b).) In appendix A the param- 
eter R(0) is found in terms of (3, and 5 is related to 0 via the condition 
of inextensibility. The resulting equations are: 


R(e) = 



and 


£5 = ± 1 X /ap\ 2 

d0 cos (3 V \<3.9 / 


(30) 


(3D 


Application of Newtonian Impact Theory 
The entire analysis up to this point is based upon equilibrium considera- 
tion of the sail subjected to normal pressure, and the equations derived are not 
dependent upon any specific aerodynamic theory. In this section use is made of 
Newtonian impact theory, (see, for example, jj?,loJ.) This theory, often 
used for hypersonic velocities, has the advantage of yielding pressures X 
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which are functions of 0 alone and thus permitting satisfaction of the boundary 
conditions (27) and (28). In addition, this aerodynamic theory has also been 
applied in p7j to a rigid idealization of a paraglider wing; therefore, a direct 
comparison can be made with the results of [7] to evaluate the use of rigid 
models in the study of paraglider behavior. 

In Newtonian impact theory, the pressure at a point is given by the 
relation 


or 


X = 2q sin 2 e if e ^ 0 


X = 0 if e < 0 






(52) 


where q is the dynamic pressure and e is the angle between the local stream- 
wise unit tangent vector and the free-stream velocity vector. Newtonian impact 
theory requires that the moving stream gives up its "normal" component of 
momentum (to the surface impacted) but retains the tangential component which 
passes off tangentially to the local surface. Only portions of the surface that 
"see" the flow have a nonzero pressure coefficient, as indicated in equa- 
tions (52). The unit normal v to the surface has been determined in terms of 
the parameters (3 and 5 in appendix A (see eq. (52)). If it is noted that 
sin e = cos(i,v) = i • v, and use is made of equation (52), the pressure on the 
sail at a given point (eq. (52)) is readily expressed as 


-.2 


X = 2q 


(cos p sin a - 


sin p cos 5 cos a) ^ 1 


f m 2 + 

,dej 


sin 6 cos a 


dp 

d9 


( 55 ) 


From equation ( 55) , it is evident that X is a function of 0 alone. 
For this case, the boundary conditions (27) and (28) on the stress resultants 



at the trailing edge of the sail are satisfied by 


(RX) = (34) 

X/jP 

Here the constant of integration C is to be determined by the relative posi- 
tions of the keel and leading edge. 


Equation (3*0 upon substitution from equation (33) for X and equation (30) 
for R yields a differential equation from which the deflected shape of the 
sail is determined; thus 



One other relation is needed. This is provided by integration of equation (3l) • 



The negative sign in equation ( 31 ) is ruled out since the lower (or upstream) 
surface must "see" the flow. Solution of the simultaneous nonlinear equa- 
tions (35) and (36), for (3 and 5, and satisfaction of specified conditions 
^( q l) an( ^ i ea( ^i n S edge of the sail, determines the deflected 

shape of the sail. This solution is effected by using finite differences to 
represent derivatives of p with respect to 0: 
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where A is the spacing of the stations and station n = 0 is taken at the 
keel. Parabolic integration is employed to evaluate the integral in equa- 
tion (36) so that 


^n ~ ^n-1 + ^(5fn + 8f n _i " f n-2) 


(38) 


where 


f n = 


cos PnJ \dey n 


Then, if values of (d|3/d0)Q, (^C /ql and angle of attack a are speci- 

fied, equations (35) and (38) became two simultaneous equations for determining 
5 n and Pn+i in terms of their values at preceding stations. For each set of 

selected values of (dp/d0)Q, and a, these equations are applied 

successively to obtain a set of values of and 5(9^) at the leading 

edge. The values of (d3/d0) o and (c jql are then varied until the 
deflected shape with the desired boundary values of P(®l) an< i '“' s 

obtained. 

The stress resultants may now be evaluated. With the relations (3I) for 
RX, and (l) for x,p equations (10), (25); and ( 26 ) become 


Nq = ( — \ ^(sin 0 + A cos 0)^ 

lllK 3 )* 3 

N x » sin 6 + A cos 0)(cos 0 - A sin 0)‘ 


) ( 39 ) 


N. 


x0 - 


\^(sin 0 + A cos 0) 2 (cos 0 - A sin 0) 


V^K 


3 L3 


J 


where A is given by equation (2) and where 


C 


has been determined in the 


/ 

course of the numerical integration described above. 

Determination of the two principal stress resultants yields the value zero 

and 


N = JL_ £i(l + A 2 )(sin 6 + A cos 0) (40) 

ol y? A 

acting perpendicular and parallel to the trailing edge, respectively. Prom fig- 
ure 2(a) and equation (l) it can he seen that lines parallel to the trailing 
edge are given by 


( X ^9r=Q A 

sin 0 + A cos 9 


(41) 


so that along these lines the maximum principal stress resultant is a constant: 


N = 


(x) 


~^(l + a2 ) 


(42) 


Integration of the stress resultants given by equations (39) at the keel 
and leading edge yields the resultant forces applied by the sail to the booms. 
(See fig. 4.) Then by considering the components of the boom forces in the Zq 
and Xq directions, the lift and drag forces are obtained. The boom forces 
and the lift and drag forces are determined in appendix B. 


NUMERICAL RESULTS AND DISCUSSION 

Numerical results are presented and discussed for a configuration that 
corresponds in all possible respects to a rigid idealization analyzed 
in |^7j • This configuration has keel and leading-edge booms of equal length; 
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and, in the undeflected flat condition, the angle 0 l (between the keel and 
leading-edge booms) is 45°. Boundary conditions which locate the leading-edge 
boom relative to the keel boom are chosen to correspond to the rigid idealiza- 
tion identified in [ 7 J as having "l80° canopy inflation." These boundary 
conditions are given by 


e (®0 =0 

5(e L ) = 28.2° 

which yields a distance between the ends of the leading-edge and keel booms 
equal to 

d = Z K ^ 2-2 cos (3^0 l ^cos 8^0 l ^ 

= 0.4872Z k 

The differential equation (35 ) , for the deflected shape, and its auxiliary 
equation ( 3 6 ) were solved with finite differences applied at intervals of 
A * 1° in the variable 0. The problem was programed on a high-speed digital 
computer, and deflected shapes were calculated for various angles of attack. 
From the calculated deflected shapes, the pressure coefficient, the stress 
resultants, the resultant forces applied by the sail to the booms, the lift and 
drag coefficients, and the ratio of lift to drag were calculated for each angle 
of attack. These results are given in table 1 and figures 5 through 7- The 
deflected shape of the surface of the sail is shown by the solid curves in fig- 
ure 5 for angles of attack of 30°, 60°, and 90° • Comparison with the rigid 
idealization of £ 7 ] (the dashed curve) shows that the rigid idealization 
is substantially in error at all angles of attack. The variation of pressure 
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coefficient Cp = X/q with angle of attack is shown in figure 6. Figures 5 
and 6 show that, as the angle of attack increases, the point of maximum deflec- 
tion and the point of maximum pressure coefficient move outboard from the keel. 

At any point on the sail, the directions of the principal stress resultants 
are perpendicular and parallel to the trailing edge. The first one has the 
value zero for all points on the sail. The other has a constant value N for 
all points on a line parallel to the trailing edge. This value is obtained by 
use of equation (42) 


N ' ^ 7 W x) ^ 

where (x)^^ gives the intersection of the line with the keel. The magnitude 

n 

varies with angle of attack through the constant which can be obtained by 

use of table 1. The magnitude N is also proportional to the distance (x) 

u 

from the nose. Note that the maximum stress is at the trailing edge and is a 
quantity of principal interest in the design of the sail. 

The lift and drag characteristics of the flexible wing are compared with 
those for the rigid idealization in figures j(a ) and 7(b) which show that the 
arbitrary shape used for the rigid idealization causes an overestimation of 
the lift and drag forces over the complete range of angle of attack; however, 
the same trends are apparent for both the rigid and flexible wings. The lift- 
to-drag ratio (fig. 7(c)) which is a measure of the angle of glide, shows a very 
significant difference for angles of attack below 50°. Computations were not 
carried below 25° because of the large amount of computer time needed in order 
to obtain the deflected shape in this range of angle of attack. This difficulty 
arises because at low angles of attack small changes in the assumed values of 
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C jql ^ and (d P /de) 0 produce large changes in the deflected shape. Presumably, 
L/D drops abruptly at some lower value of a. But the results already show 
that L/D for the flexible sail is more than three times the maximum value 
of L/D for the rigid idealization. 

The resultant forces, applied by the sail to the keel and leading-edge 
booms, and the location of these forces are given in table 1. These values may 
be used by the designer for calculating stresses in shroud lines and spreader 
bars. Stresses in the booms may then be found by analyzing these booms as 
free-free beams subjected to the self-equilibrating system of loads provided by 
the sail, shroud lines, and spreader bars. 

Effects of changing the dihedral of the wing were also considered. The 
boundary conditions for this case correspond to maintaining the value 

d = 0.4872 Z k 

and specifying various values of These computations were performed for 

an angle of attack a = 35° in order to find deflected shapes, pressure dis- 
tributions, boom forces, and lift and drag coefficients. The variation of the 
force characteristics with dihedral is given in table 2. The variation of lift 
and drag coefficients with dihedral is shown in figure 8. The maximum values 
of lift and drag coefficients are obtained for a dihedral angle of about -8°. 

CONCLUDING REMARKS 

Equilibrium equations for the sail of the paraglider wing have been derived 
and integrated. Results are first obtained in rather general form for the 
stress resultants in the sail, boom forces and lift and drag in terms of the 
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pressure on the sail and a parameter that describes the deflected shape of the 
sail. The specification of an appropriate aerodynamic theory - in the present 
case, Newtonian impact theory - then permits satisfaction of the boundary condi- 
tions at the trailing edge and calculation of the deflected shape of the sail. 
Finally, the general formulas are applied to calculate the stress resultants, 
lift, and drag. 

Numerical results have been obtained and have been compared with published 
results for a rigid idealization of the paraglider wing. The comparison shows 
that the assumed shape of the rigid wing is considerably in error over the 
complete range of angles of attack. Consequently, the lift and drag coeffi- 
cients and, especially the lift-to-drag ratio, for the flexible wing are sig- 
nificantly different from the values found for the rigid wing. Thus use of 
rigid idealizations in wind-tunnel investigations to draw conclusions as to the 
aerodynamic characteristics of the paraglider - especially the lift-to-drag 
ratio - may be greatly misleading. 

The calculated stress resultants and boom forces provide a basis for 
design of sails, booms, shroud lines, and spreader bars for a paraglider in 


hypersonic flight. 
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APPENDIX A 

SHAPE OF LOADED SURFACE OF THE SAIL 


In this appendix the shape of the loaded surface of the sail is represented 
by the angles (3 and 6 shown in figure 2(b)). The x-j_-axis of the rectangular 
Cartesian coordinates x-^ y , z^ is alined with the keel of the paraglider 
and the x-^z-^-plane is the vertical plane of symmetry of the paraglider. Hence 

z^_ = x sin 3 

x^ = x cos 3 cos 5 > ( 43) 

y-j_ = x cos 3 sin 5 j 

Consider the keel at an angle of attack a and define rectangular 
Cartesian coordinates x Q , y Q , z Q such that the XQ-axis is alined with the 
direction of airflow. The XQZQ-plane is coincident with the x^Z]_-plane (fig. l) , 
thus 


Zq = x(sin 3 cos a - cos 3 cos 5 sin a) 

Xq = x(sin 3 sin a. + cos 3 cos 5 cos a) > (44) 

y Q = x cos 3 sin 6 j 


Now if zq is treated as a function of Xq, y^, representing the deflected 
surface in the Xq, y^, Zq coordinate system. 


Sz n Sz n 

dz Q = — 0 dx Q + -0 dy Q 
oxq oy-f 


(^5) 


'0 


and it can be shown from equations (44) and (45) that: 
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> 7 -sin 6 cos a ^ + (sin B cos 6 cos a - cos B sin a) cos 6 — 
oz, 0 _ d9 _ d9 

^*0 -sin 5 sin a ^ + (sin 3 cos 5 sin a + cos 3 cos a) cos 3 — 

d0 d0 


(46) 


and 


^0 


cos 5 ^ + sin 3 cos 3 sin 5 — 

§ de (4?) 

-sin & sin a — + (sin 3 cos 6 sin a + cos 3 cos a) cos 3 — 
d0 d0 


d5 


Also the square of the length of a line element on the deflected surface is 
given by 


ds 2 = cLxq^ + + dz, 0 C 


(48) 


which upon substitution from equations (44) and (45) yields 

,2 . / acsN 2 | 


ds 2 = dx 2 + x 2 


(I) + ^(i) 


de £ 


(49) 


But inextensibility requires that the metric (eq. (5)) remain unchanged; thus, 
it follows from equation (49) that: 


= ± i /q _ /d£ 

d9 cos 3 V \d9 


(50) 


The unit outward normal to the deflected surface is given by 


v = ± 


t ft of ♦ 1 


OX, 


0, 


dy r 


1/2 


8z n - 8z n _ 

— 2 . i - _ a j + k 

oxq oy 0 


(5D 


where the positive sign refers to surfaces which are concave downward, and the 
minus sign to those which are concave upward. The unit normals i, and k 

are directed along the x Q -, y Q -, and z Q -axes as shown in figure 1. Upon sub- 
stitution from equations (46), (47), and (50) for ^2., ^2., and — 

6 xq 6y 0 *9 


and 
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equation (51) yields 


v = 


(cos 0 sin a - sin fl cos 6 cos a) i/l - + sin 8 cos a 

v y \de / de 


sin p sin Si/l - ( + cos 5 ^ 

d9 ) de 


J 


(cos 0 cos a + sin 0 cos 5 sin a),/l - \ - sin 5 sin a — 

V \ de/ de 


(52) 


The unit tangential vector in the x-direction is given hy 


t-j^ = (sin p sin a + cos p cos 5 cos a)i 

+ (cos P sin 8)j + (sin P cos a - cos p cos 5 sin a)k 


(53) 


and the unit tangential vector in the 6-direction hy 


Tg = v x t-l 


-sin 8 cos 


cos 8 


a ./l - + (cos p 

v \dey 

F( 


sin a - sin p cos 8 cos a) — 

de 


— ^ - sin p sin 8 — 

de/ de 


sin 8 sin a Jl - ( ^ 

de 


^ + (cos p cos a + sin p cos 8 sin a) — 

/ de 


(5*0 


The quantity g can be expressed in terms of v and Tg [ll] and then in 
terms of the angle p by using equations (52) and (5^): 


g = XT o * v 


2 v ,2 


Jl - (&) 

» \cLe / 


I sin P 

2 V 08 3 


-(H) £ 


de^ 


(55) 
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Then from equation (6) the function R(9) can be expressed in terms of the 
quantity 3 as follows: 
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APPENDIX B 

CALCULATION OF LIFT AND DRAG FORCES ON THE PARAGLIDER WING 


In this section the lift and drag forces are derived. As a preliminary 
step, the resultant forces applied hy the sail to the keel and leading-edge 
booms are obtained. These resultant forces and their locations are useful in 
obtaining the forces in the shroud lines and spreader bars of the paraglider. 

The resultant force applied to the keel by the half of the sail considered 
in the analysis is obtained by integration of the stress resultants for the 
sail along the keel boundary. Since the tangent vectors along the keel do not 
vary with x the vector equation for the force is given by 



and this force acts through the point x = x^- of the keel, where 


X K 



(58) 


Here and {^ 2 } are tangent vectors of the surface at the 

keel (eqs . ( 53 ) and (54) evaluated at 0=0). Similarly, the vector equation 
for resultant force applied by the sail to the leading-edge boom is given by 



and this force acts through the point x = of the leading edge where 
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X L = 


JO 


; ("e)e=e L dx 


(60) 


f dx 

0 ' ej 0=0L 

In equation (59) (^l) anc ^ (^ 2 ) are un ^ tangent vectors of the surface 

(eqs. ( 55 ) and ( 5 * 0 ) at the leading edge 0 = 0 ^. 

It is necessary now to determine the components of f k and F^ in the 
x G , y Q , and z 0 directions. For example, the component of F^ in the 
x Q direction is given by 


F Kx 0 = P K 


1 = cos a 


'•K 


( NxQ )o^) dx + s ^" n a 


M\ f h („ e ) 

de 4=o J o ' ' e=<J 


dx 


(61) 


or for Newtonian impact theory (eq. (39)) in nondimen si onal form: 


Fl&o _ ! 

i C ^ 

l K) 

ll L 

qs 2 1 


sin 0 ^ 


— cos 
A 


a - ( sin a 

We, 


'9=0 


(62) 


Similarly the remaining components are 


* K yo _ 1/ C \ V' 


qS 


2 qZ K 3 Sin 0 L 


- l^-f 

\ dQ /Q=0 


(63) 


fK2 o _ 1 / c \ z k/ z l 

qS 2| 


Wk 


3 sin 0 L 


— sin 


a + ( cos a 

\ d0 ) e =o 


(64) 


*LXq 

IS" 


= i /c \( 1 k/ 1 l ) 2 / t 
' 2 ^k 5 ) sln 0L \ 1 


cos e L - 
sin ©l 


- ^ 2 


(65) 
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where 


■Ly 0 = 

qS 


/ c \ 

.(W 1 l) 2 | 

( 

J-i 

cos 0 L - l L /l K 

(i'k 3 / 

) sin 0 L 1 


sin 0^ 

/ c ' 

( 1 k/ 1 l) 2 | 

Kn 

cos 0 L - T'i J lT'fc 

\ ql K 5 y 

1 sin 0 L 

V 

sin 0^ 


- V 


Ij_ = sin p^Qj-^sin a + cos p^e^^cos S^ 0 L ^cos a 
J-L = cos P^Q-j^sin 6 ^ 0 ]-^ 

K-^ = sin p^Qj^cos a - cos p^e L ^cos S^ 0 L ^sin a 


12 = -sin 5^0-^cos a ^1 - { — ^ 

+ jcos P^Qj^sin a - sin p^L^cos 5^0^ 


d0 /0=0 L 


cos a 


dp\ 

de /0=i 


0=0 L 


J 2 = cos 6 


nRu ■ sln ^ L)sia 6(0 l) (^L l 


(m e 


K 2 = sin 8 (e L )sin a 1 - (SE 


\ d0 ^0=0T 


|cos P^l^cos a + sin p^B^^cos 5^0^ sin a 


d %- 


0 ~ 0T 


The coordinates Xq^, etc. (fig. k) through which these forces act 
written with the aid of equations (44 ) , ( 58 ), and ( 60 ) as 


(66) 

(67) 

( 68 ) 

( 69 ) 

( 70 ) 

(71) 

(72) 


(73) 
may be 
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X 0K “ 3 l K cos a 


y 0K = ° 


Z 0K = - 3 l K sin a 


X 0L " 3 Z L I 1 


2 

'OL 3 


y^r - - z l j i 


( 7*0 

( 75 ) 

( 76 ) 

( 77 ) 


( 78 ) 


2 

z OL = 7 Z L K 1 


( 79 ) 


Finally, when both halves of the symmetrical sail are considered (fig. 4) , the 
lift and drag coefficients can be expressed as follows: 




+ F 


LZr 


(80) 


C 


D 


D 

qS 


qS ( F K x 0 + F Lxq) 

where L and D are the lift and drag forces, respectively, on the wing; 
1 P 

q = -pV is the dynamic pressure, and 


S = ZprZ^ sin 9 ^ 


(82) 
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is the total surface area of the sail. The resultant of the lift and drag 
forces acts through the point (xq, where 


x 0 = 


x OL F Lz 0 + x OK F Kz 0 
F Lzo + F Kz q 



Equations ( 80 ) and (8l) have been used to obtain the results shown in figures 7 
and 8. 
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TABLE 1.- FORCE CHARACTERISTICS OF A PARAGLIDER WING 


AT VARIOUS ANGLES OF ATTACK, a 



— = i; 6 
_ L 

» 4-5°; P( 

'e L ) - 0; 

6(6 L ) = 28.2°; iL 

= 0.487; 

a, deg 

25 

50 

55 

40 

45 

50 

(dp/d0) o 

0.99895 

0.99605 

0 . 99254 

0.98854 

0.98455 

0.9804 

c/ql K 5 

. 004155 

.01046 

.01979 

.05218 

.04716 

. o645< 

FKxq/hS 

. 000157 

. 00105 

.00522 

. 00724 

.0154 

.0221 

F Ky 0 / qS 

.000154 

.000658 

. 00175 

.00546 

.00588 

.0089; 

%z 0 / qS 

. 00516 

.00791 

.0147 

.0255 

•0550 

.0452 

p Lxo/ qS 

. 000791 

. 00244 

.00548 

.0105 

•0175 

.0265 

F Ly 0 / qS 

- .00244 

- . 00578 

-.0105 

-.0159 

-.0221 

- . 0289 

F i-o/ qS 

. 00185 

.00497 

.00967 

■ 0157 

.0227 

.0500 

x ok/ z k 

.6o4 

•577 

. 546 

.511 

.471 

.428 
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TABLE 2.- FORCE CHARACTERISTICS OF A PARAGLIDER WING 
AT VARIOUS DIHEDRAL ANGLES, p(e L \ 


^ = 1; 0 L = 45°; A = 0.4872; a = 35° 
_L 


p(e L ), deg 

-15 

-10 

-5 

0 

5 

10 

i 4.4 

(dp/d0) Q 

0 . 98727 

0.98702 

0.98905 

0.99234 

0.99574 

0.99842 

0.99970 

c/qI K 5 

. 01894 

.02126 

•02153 

•01979 

.01645 

.01190 

.007500 

W qs 

.00304 

.00341 

.00347 

.00322 

.00270 

.00196 

.00124 

F Ky 0 / qS 

. 00213 

. 00241 

. 00225 

.00173 

.00107 

.000473 

.000130 

F Kz 0 / qS 

.0140 

•0157 

.0160 

.0147 

.0123 

.00888 

.00560 

W qs 

.00778 

.00787 

.00701 

.00548 

.00366 

.00192 

.000760 

v ./* 5 

- . OO633 

- . 00859 

-.0100 

-.0103 

-.00935 

-.00723 

- . 00475 

P LZ 0 / qS 

.0105 

.0114 

.0110 

.00967 

. 00760 

.00520 

. 00313 

x ok/ z k 

.546 

.546 

• 546 

. 546 

.546 

.546 

.546 

f ok/ 1 k 

0 

0 

0 

0 

0 

0 

0 

z ok/ 1 k 

-.382 

-.382 

-.382 

-.382 

-.382 

-.382 

- • 382 

x ol/ 1 k 

.383 

.415 

.448 

.481 

•515 

.548 

• .576 

^ o ~ lI 1 k 

.264 

• 293 

• 310 

• 315 

.310 

• 293 

.368 

z ol/ 1 k 

-.478 

-.431 

-.384 

-•337 

-.289 

-.242 

-.201 

C L 

.0490 

.0542 

.0540 

.0487 

.0397 

.0281 

•0175 

C D 

.0216 

.0226 

. 0210 

.0174 

.0127 

.00777 

.00400 

l/d 

2.26 

2. 4 o 

2.57 

2.80 

3.12 

3.62 

4-37 

x oj l K 

.476 

.491 

.506 

.520 

• 534 

.547 

•557 

z o/*K 

-.451 

-.417 

-.384 

-.354 

-.329 

-.313 

-.314 





Figure 1.- Paraglider wing configuration. 








Figure 2.- Concluded- 
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Figure 7*- Concluded. 
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Figure 8.- Variation of lift and drag characteristics of paraglider 
wing with dihedral angle |3(e L ). Angle of attack, a = 35 • 
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